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摘 ” 要, 基于 子午 工程 项 目 北京 钠 激 光 雷 达 的 长 期 观测 数据 (2010 年 4 月 到 2011 年 9 月 ), 采用 
Yang 等 人 (2008) 的 重力 波 参 数 提取 方法 ,提取 了 在 253 个 夜晚 2 208 h 的 有 效 观 测 时 间 里 出 现 的 162 
个 钠 层 单 色 重力 波 的 垂直 波长 .周期 、 幅 度 值 和 幅度 增长 因子 参数 值 。 统 计 分 析 显 示 激 光 雷 达 所 观测 
到 的 中 层 顶 区 域 重力 波 县 有 系统 性 的 参数 关系 ,A.=0.2267T%”” .KE=6.66x10-* 有 和 KE=8.97x10* x 
凡 ”“, 这 些 实验 观测 结果 验证 了 Gardner(1994) 扩 散 滤 波 理论 的 结果 ,但 同时 发 现 仅 认为 波 致 涡流 扩散 
是 重力 波 衰减 的 主要 途径 是 不 完全 正确 的 ， 其 它 因素 引起 的 衰减 在 重力 波 的 饱和 和 耗 散 机 制 中 也 起 
到 了 重要 作用 。 
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Abstract: Based on the long -term observations from Beijing sodium lidar of Chinese Meridian Project, 
162 monochromatic gravity waves was distinguished in the sodium layer during the 2208 hours’ 
observation in 253 nights (from 2010 April to 2011 September). The vertical wavelength, wave period, 
wave amplitude and its growth factor of every monochromatic gravity wave were precisely extracted with 
the method derived by Yang et al.(2008). Statistical results reveal that there are systematic relationships 
between those wave parameters, and they respectively are A.=0.2267T%;”, KE=6.66x10™™" ka2"™ and KE= 
8.97x10™ fis"™. Those results partially agree well with the results predicted by Gardner ’s(1994) diffusive 
filtering theory. However, it is found incorrect to regard the wave—induced eddy diffusivity as the only 
main damping processes of gravity waves. It is suggested that the damping processes caused by other 
possible factors also play important roles in the saturation and dissipation mechanisms of gravity waves. 
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0 引 塞 


大 气 重力 波 活动 是 一 种 大 气 中 常见 的 物理 现 
象 , 大气 的 温度 、 密 度 、 风 场 等 特征 参量 会 随 着 重力 
波 活动 发 生 振荡 起 伏 ; 同时 重力 波 活动 将 能 量 和 动 
量 从 低层 大 气 传输 到 中 间 层 上 部 及 低热 层 中 ， 对 大 


气 中 大 尺度 物质 、 能 量 和 动量 的 交换 和 重新 分 配 有 
着 重大 影响 ,中 层 顶 区 域 重力 波 活动 频繁 ,通过 重力 


波 的 传播 实现 低层 大 气 和 高 层 大 气 的 耦合 ; 对 此 区 
域 中 重力 波 活 动 的 研究 和 籽 有 利于 对 大 气 中 的 物质 、 
能 量 和 动量 的 传播 、 沉 积 以 及 收 支 平衡 的 动力 学 过 
程 的 理解 由， 

由 于 钠 层 中 的 Na 原子 是 中 层 顶 区 域 重力 波 活 
动 很 好 的 示 踪 物质 ,有 上 县 有 高 时 间 和 空间 分 辩 率 的 Na 
激光 雷达 是 研究 该 区 域 重力 波 的 传播 、 饱 和 和 和 耗 散 
特征 的 有 效 工 具 中。 基于 激光 雷达 的 长 期 观测 ， 
Gardner(1994) 等 研究 了 重力 波 的 参数 关系 ,并 提出 
了 扩散 滤波 理论 外。 该 理论 认为 波 臻 涡流 扩散 效应 
是 重力 波 耗 散 的 主要 途径 ， 重 力 波 在 向 上 传播 的 过 
程 中 会 被 背景 大 气 选择 性 吸收 ， 只 有 那些 涡流 扩散 
速度 小 于 垂直 相 速 度 的 波 才 能 向 上 传播 ， 重 力 波 垂 

直 波 数 谱 相 应 地 出 现 选择 性 滤波 效应 ; 并 且 由 此 得 
出 雷达 观测 到 的 中 层 顶 区 域 重 力 波 是 趋 近 饱和 的 ， 
县 有 系统 性 参数 关系 A.xT%” KExk” 和 KE cf 。 
然而 ， 重 力 波 的 线性 饱和 理论 认为 重力 波 垂直 波 数 
谱 的 形状 和 幅度 值 是 由 风 剪 切 和 对 流 不 稳定 性 控制 
的 ， 饱 和 重力 波 的 垂直 波 数 谱 在 高 波 数 区 域 的 谱 指 
数 近 似 为 n=-3!H。Hines(1991) 的 多 普 勒 扩展 理论 和 
Weinstock 的 尺度 相关 的 扩散 理论 则 分 别 认 为 重力 
波 的 饱和 和 耗 散 完全 由 波 - 波 相互 作用 的 非 线 性 效 
应 和 混沌 扩散 运动 控制 "9 。 

基于 国家 “子午 工程 项 目 ” 的 北京 激光 雷达 站 的 
大 量 观测 数据 ， 文 章 采 用 Yang 等 的 重力 波 参 数 提取 
方法 外 ,提取 了 162 个 Na 激光 雷达 观测 到 的 单 色 重 力 
波 的 垂直 波长 周期 .幅度 值 和 幅度 增长 因子 等 参数 
值 ,分 析 了 北京 上 空 雷达 观测 到 的 中 层 顶 区 域 重力 波 
的 参数 关系 ,从 而 验证 了 大 气 重力 波 的 相关 理论 。 


1 钠 层 重力 波 参 数 提取 方法 


基于 大 气动 力学 基本 原理 ,Gardner and Voelz 给 


出 了 钠 层 重力 波 扰动 和 参数 之 间 的 关系 叫 . 
no{z—yHIn[1+[Ae®/(y—1)]cos(wt—k:7)]} (1) 
{1+[Ae®/(y—1)]cos(wt—k: 7)]} 
式 中 :n,(7,) 为 钠 层 密度 对 重力 波 的 响应 ;n 为 本 底 
钠 层 ;Aes 为 重力 波 在 高 度 z 处 的 幅度 值 ;8 为 幅度 
增长 因子 ;w 为 重力 波 频 率 ; k=kXx+k.2 为 波 数 矢量 ，; 
万 为 大 气 标高 ;y 为 热 熔 比 ; 置 坐标 ;t 
为 时 间 。 
在 重力 波 的 幅度 不 太 大 的 情况 下 ， 仅 考虑 重力 
波 的 线性 扰动 , 式 (1) 可 以 转换 为 式 (2)， 


ns(r,D) ~no(z)-[noz)+YH a 2 = 4 


由 于 重力 波 引起 的 钠 密 人 
直接 相关 ,相对 钠 密度 扰动 (rsam) 可 以 用 大 气 密 度 扰 
动 (rammee) 表 示 为 叫 : 


71,(7, 门 = 


r=xX+ZZ 为 位 


cos(wt—k:r) (2) 


1+78 dm 


no de ocos(wt kr) (3) 
y—1 


采用 Yang (2008) 等 的 方法 ， 利 用 Levenberg-- 
Marquardt 拟 合 方法 ， 找 到 较为 真实 的 本 底 钠 层 m， 
计算 出 重力 波 活动 导致 的 钠 层 密度 扰动 An,=(n,-n0)/ 
m, 就 能 由 式 (3) 得 到 重力 波 线性 扰动 对 应 的 大 气 密 


度 扰动 一 二 Aescos(wt-k.r); 通 过 测量 大 气 密度 拓 
ee 就 可 以 得 到 重 


力 波 的 垂直 波长 A.; 同时 对 重力 波 扰 动 廓 线 的 波 腹 
作 线 性 拟 合 log(Z)=log(4)+6.z ,就 可 以 得 到 重力 波 
在 不 同 高 度 上 的 幅度 值 4e#, 也 可 得 到 重力 波 在 传播 
中 的 幅度 增长 因子 6%。 另 外 ,在 密度 廓 线 序列 中 测 
量 波 峰 ( 谷 ) 向 下 的 移动 速度 ， ee 
直 相 速度 c.。 重 力 波 的 观测 周期 Tw 就 可 以 通 

和/T 计算 出 来 外 。 

对 于 大 多 数 低频 重力 波 而 言 ,在 满足 |2kH|>>1 
的 情况 下 ,由 重力 波 的 非 线 性 色散 关系 , 单 色 重力 波 
KE 动能 可 以 计算 为 中 . 


KE(K， 中 )= | 于 3 lw = ey (4) 


其 中 v. 和 v 分别 为 a 
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直 和 水 平方 向 上 的 风速 。 

需要 说 明 的 是 ， 使 用 上 述 方法 准确 提取 重力 波 
参数 时 ,需要 注意 下 列 四 个 前 提 假 设 乌 :(]1) 钠 层 密 度 
分 布 具有 水 平均 匀 性 ;(2) 有 重力 波 活动 时 , 钠 层 化 
学 反应 对 钠 层 的 影响 很 小 ;(3) 忽 略 重 力 波 扰动 引起 
的 扩散 和 化 学 反应 效应 变化 ;(4) 除 了 化 学 反应 和 钠 
层 水 平 运动 外 ， 钠 层 的 变化 主要 是 由 重力 波 扰动 引 
起 的 .其 中 ,(1)-(2) 假 设 主 要 是 用 来 排除 钠 层 中 出 现 
ee 这 在 重力 波 参 
数 提取 时 ， 通 过 观察 钠 层 形态 变化 就 能 做 出 正确 判 
断 。 ee ,Xu 的 理论 模型 研究 都 指出 ， 
由 于 化 学 反应 的 时 间 周 期 较 长 ， 在 研究 重力 波 扰动 
时 , 钠 层 化 学 反应 效应 变化 的 影响 可 以 不 考虑 中, 实 
际 情况 下 ， 化 学 反应 和 水 平 不 均匀 性 的 影响 是 必然 
存在 的 ;所 以 ,在 重力 波 参数 提取 时 ,需要 根据 钠 层 


形态 变化 的 实际 情况 正确 估计 这 些 因素 对 重力 波 参 
数 准确 性 的 影响 。 
2 激光 雷达 观测 结果 


从 2010 年 4 月 到 2011 年 9 月 ,在 北京 Na 激光 
雷达 的 253 个 夜晚 (2208 hb) 的 有 效 观 测 数据 里 ,共计 
观测 到 了 162 个 单 色 重力 波 事件 。 基 于 上 述 的 重力 
波 参 数 提取 方法 ， 提 取 了 每 一 个 重力 波 的 垂直 波长 
A.、 周 期 7 .幅度 增长 因子 6 和 幅度 值 4es% ,并 在 
1 中 对 应 地 画 出 了 出 现 的 时 间 分 布 图 。 观 测 到 的 
重力 波 波长 的 平均 值 为 3.73 km， 周 期 的 平均 值 为 
3.38h ,在 90km 处 的 平均 幅度 值 为 2.94%。 


N=162 Mean=3.73 km | 
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1 激光 雷达 观测 到 中 层 顶 区 域 重力 波 的 参数 值 及 对 应 出 现时 
间 分 布 图 


Fig.1 Distribution of the wave parameters versus the day number of 


哆 


year, which is observed in the mesopause region by lidar 


2 的 统计 显示 ,重力 波 的 垂直 波长 和 A. 分布 在 2~ 
11.88km ,常见 的 重力 波 波 长 为 2~4km; 对 应 的 周期 T 
分 布 在 1.03 ~14.97h ,常见 的 周期 为 1~4h; 重 力 波 的 幅 


度 值 都 小 于 20% ,常见 的 重力 波幅 度 值 小 于 4% 。 
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2 激光 雷达 观测 到 的 重力 波 的 垂 
和 幅度 值 的 统计 直方 图 


Fig.2 Histograms of the vertical wavelength , observed period , growth 
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factor ,and wave amplitude for AGWs observed by lidar 


基于 上 面 得 到 的 各 个 重力 波 参数 值 ， 对 北京 上 
空 激 光 雷 达观 测 到 的 中 层 顶 重力 波 的 参数 关系 作 了 
如 下 的 详细 分 析 。 

3 是 重力 波 垂直 波长 A 与 周期 Tw 之 间 的 对 数 
坐标 关系 ,可 以 看 出 两 者 具有 很 强 的 系统 相关 性 。 对 
图 中 这 些 数 据点 作 .=C7% 的 最 大 似 然 拟 合 ， 发 现在 
北京 观测 到 的 重力 波 有 7% 吕 对 垂直 波长 A. 的 系统 相 
关 性 (虚线 ), 即 A=0.2267%™” ;尽管 认为 背景 大 气 风 场 的 
Doppler 效应 会 给 周期 Tv 的 测量 带 来 较 大 的 Doppler 
偏差 ,其 拟 合 相关 系数 仍 达 到 了 六 =0.819。 同 时 ,分 别 
对 出 现在 夏季 (3~8 月 份 ) 和 冬季 (9~12 月 ,1~2 月 ) 的 重力 
波 作 了 统计 分 析 ; 在 两 个 季节 里 ,垂直 波长 和 .与 周期 7 
之 间 的 关系 分 别 是 .=0.2457%" 和 A.=0.2017%* 。 从 
夏季 到 冬季 ,观测 到 的 Cp 值 变化 较 大 ,这 说 明 在 北京 
地 区 重力 波 活 动 的 季节 变化 较 明 显 。 
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3 激光 雷达 观测 到 的 重力 波 垂直 波长 长 A. 与 周期 7 之 间 的 
系统 关系 。 虚 线 是 对 数据 点 的 A.=C7T% 最 大 似 然 拟 合 ,其 中 
C 是 拟 合 得 到 的 常数 系数 ,p 是 斜率 。 记 是 拟 合 相关 系数 


Fig.3 Systematical relations between the vertical wavelength A, and 


period T,, of gravity waves observed by lidar. Dash lines are 
the maximum likelihood fits to the data in form of A.=CT,,, 
where C is the constant coefficient and p is the slope of 


logarithmic plot. ris the correlation coefficient of fitting 


4 和 图 5 分 别 给 出 了 观测 到 的 重力 波 在 90 km 
处 的 KE 动能 与 波 数 (k/27) 和 频率 (f,) 之 间 的 关系 ， 
发 现 KE 动能 和 波 数 之 间 有 KE=6.66x10-k* 中 的 系 
统 相 关 性 ,KE 动能 和 频率 之 间 有 的 相关 性 KE= 
8.97x10 人 a” , 且 这 些 参 数 关系 在 夏季 和 冬季 也 略 有 
不 同 。 
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4 激光 雷达 观测 到 重力 波 的 KE 动能 与 波 数 Coy/2m 的 系统 关系 。 
虚线 是 对 数据 点 的 .=C7% 最 大 似 然 拟 合 ,其 中 C 是 拟 合 得 
到 的 常数 系数 ,p 是 斜率 。 产 是 拟 合 相关 系数 

Fig.4 Systematical relations between the kinetic energy and wave 
number of gravity waves observed by lidar. Dash lines are the 
maximum likelihood fits to the data in form of A.=CT,, , where 
C is the constant coefficient and is p the slope of logarithmic 


plot. Pis the correlation coefficient of fitting 
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5 激光 雷达 观测 到 重力 波 的 动能 与 频率 f 的 系统 关系 。 虚 线 
是 对 数据 点 的 A.=C7T% 最 大 似 然 拟 合 ,其 中 C 是 拟 合 得 到 的 
常数 系数 ,p 是 斜率 。 户 是 拟 合 相关 系数 


Fig.5 Systematical relations between the kinetic energy and frequency 


of gravity waves observed by lidar. Dash lines are the maximum 
likelihood fits to the data in form of A.=CT,, ,where is C the 
constant coefficient and p is the slope of logarithmic plot. rr is 


the correlation coefficient of fitting 


结果 讨论 


从 激光 雷达 的 观测 数据 里 准确 提取 重力 波 的 参 
数值 ,找到 钠 层 密度 分 布 的 真实 本 底 是 至 关 重 要 的 。 
Gardner and Voelz 在 认为 本 底 钠 层 密度 呈 高 斯 分 布 
的 基础 上 ， 采 用 计算 重力 波 垂 直 波 数 谱 的 方法 提取 
如 垂直 波长 .振幅 增长 因子 等 外 ,但 

是 ,由 于 钠 源 和 钠 汇 是 动态 变化 的 , 且 受 大 气 输 运 过 
程 影 响 ,本 底 钠 层 形状 是 不 断 变 化 的 , 且 常 常 偏离 高 
斯 分 布 的 .所 以 ,采用 Gardner and Voelz 的 方法 ,会 给 
重力 波 人 参数 的 提取 结果 带 来 巨大 偏差 由。 这 里 采用 
Yang(2008) 等 的 参数 提取 方法 ,利用 平均 钠 层 作为 
初始 本 底 钠 层 ， 再 采用 Levenberg-Marquardt 拟 合 方 
法 ,逐步 逼近 真实 本 底 钠 层 。 这 种 重力 波 参 数 提取 方 
法 是 对 Gardner and Voelz 方法 的 重大 改进 准确 
地 给 出 本 底 钠 层 和 重力 波 参数 。 

对 于 单 色 重力 波 事件 的 确定 ， 这 里 采用 Beatty 
等 所 描述 的 标准 趾 .(1) 钠 密度 廓 线 呈 明显 的 波状 扰 
动 结构 ;(2) 这 种 扰动 结 续 出 现在 3~5 个 钠 
密度 廓 线 中 , 即 出 现时 间 超 过 10min;(3) 这 些 波 状 扰 
动 结构 具有 较 一 致 的 向 下 传播 相 速 度 。 在 北京 所 观 
测 到 的 162 个 单 色 重 力 波 事件 都 符合 上 述 条 件 。 当 
然 ， 对 于 那些 近似 水 平 传播 的 重力 波 和 明显 出 现 


至 少 连 
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Doppler 频 移 的 重力 波 来 说 ,这 个 判 据 会 给 重力 波 事 
件 的 数目 统计 带 来 一 些 误差 ,并 且 ,雷达 观测 到 的 重 
力 波 周 期 7 也 应 该 是 受到 钠 层 背景 风 场 的 Doppler 
频 移 效 应 影响 后 的 周期 ， 而 不 是 重力 波 本 身 的 固有 
周期 。 但 是 ,Beatty 等 人 认为 ,由 于 雷达 观测 到 的 重 
力 波 的 波长 和 周期 往往 具有 明显 的 系统 相关 性 (如 
图 3 所 示 ) ,这 些 重 力 波 应 该 很 少 受到 钠 层 背景 风 场 
的 Doppler 频 移 效应 影响 中 ,在 完成 激光 雷达 和 
Radar 观测 结果 的 对 照 后 ,Collins 等 人 也 认为 雷达 在 
中 层 顶 区 域 所 观测 到 的 都 是 没有 明显 受到 背景 风 场 
Doppler 频 移 效应 影响 的 重力 波 ,观测 到 的 重力 波 参 
数值 可 以 代替 重力 波 的 固有 参数 值 3， 

图 2 中 的 统计 分 析 显 示 , 随 着 波长 .周期 或 幅度 
值 的 增加 重力 波 出 现 的 次 数 急剧 减少 ， 这 说 明 在 北 
京 上 空中 层 顶 附近 雷达 观测 到 的 主要 是 中 小 尺度 、 
短 周期 重力 波 。 在 相似 纬度 上 ,Illinois(40°10'N,88° 
10'W) 激 光 雷 达 小 组 也 报道 了 相同 的 观测 结果 ;在 其 
它 纬度 上 ,Yang .Beatty Collins 等 的 雷达 观测 也 给 出 
了 相同 的 结论 " 久 。 

关于 钠 激 光 雷 达 对 重力 波 的 观测 范围 ,Gardner 
and Taylor 指出 , 钠 密 度 激光 雷达 对 重力 波 的 观测 精 
度 相当 高 ， 理 论 上 能 观测 到 尺度 和 激光 雷达 的 分 辩 
率 相当 的 小 尺度 重力 波 (垂直 波长 几 百 米 .周期 几 分 
钟 )。 但 是 ,由 于 激光 器 输出 能 量 的 涨 落 和 背景 噪声 
的 影响 ， 通 常 实际 能 准确 分 辩 的 最 小 尺度 重力 波 的 
波长 大 致 在 2km 波长 左右 ,周期 30 min 左右 。 同 时， 
由 于 钠 层 本 身 的 厚度 限制 ， 激 光 雷 达 不 可 能 观测 到 
大 尺度 的 重力 波 (垂直 波长 超过 15 km ,同时 周期 大 
于 5hb) 叫 。 这 里 北京 激光 雷达 所 观测 到 的 最 大 尺度 
重力 波 的 垂直 波长 为 11.88 km ， 周 期 为 898.25 min. 
这 和 Beatty (1992 ,12.8 km) ,Gardner and Voelz(1987 ， 
17km) 和 Yang(2008 ,11.6 km) 等 的 报道 结果 相 类 似 。 
最 小 尺度 的 重力 波 的 垂直 波长 为 2km， 也 和 其 它 地 
方 的 报道 相仿 。 若 需要 观测 到 中 层 顶 附近 更 小 尺度 
的 重力 波 ,不 仅 需 要 提高 雷达 的 空间 分 辨 率 ,而 且 需 
要 有 效 地 降低 背景 和 系统 噪声 。 关 于 光子 噪声 和 激 
光 器 输出 能 量 涨 落 引起 的 重力 波 参数 的 测量 误差 ， 
Gardner and Voelz(1987 ,公式 (54-55)) 已 作 了 详细 的 
分 析 凹 。 在 通常 情况 下 ,由 于 光子 噪声 和 激光 能 量 涨 
落 的 影响 ,垂直 波长 A 的 误差 小 于 1% ,重力 波 的 观 


测 周期 7 的 偏差 则 小 于 3% 。 

由 于 大 气 密度 随 高 度 的 增加 会 明显 减 小 ， 大 多 
数 重力 波 理论 都 认为 ， 随 着 高 度 的 增加 重力 波幅 度 
值 会 不 断 增长 ,直至 饱和 。 实 际 上 ,重力 波 在 传播 过 
程 中 会 受到 背景 大 气 的 吸收 、 反 射 和 耗 散 等 作用 , 重 
力 波 的 幅度 并 不 一 定 是 随 高 度 的 增加 而 增 大 的 。 图 2 
显示 ， 几 乎 一半 的 重力 波 的 幅度 值 增长 因子 B 
(Growth Factom) 是 负 值 。 这 说 明 在 中 层 顶 附近 激光 雷 
达观 测 到 的 重力 波 的 幅度 并 不 是 随 高 度 的 增加 都 增 
大 的 。Gardner 的 扩散 滤波 理论 认为 只 有 那些 有 效 垂 
直 扩 散 速 率 (KD) 小 于 垂直 相 速 度 (w/) 的 重力 波 的 
幅度 地 会 随 高 度 的 增加 而 增 大 , 即 Dh? <w; 其 他 重 
力 波 将 会 被 严重 衰减 而 消失 名 。Yang 等 在 Sio José 
dos Campos(23*N ,46"W) 用 激光 雷达 首先 观测 到 中 
层 顶 附近 重力 波 的 这 种 衰减 现象 止 ;这 里 北京 激光 雷 
达 对 中 层 顶 附近 重力 波 活动 的 研究 也 发 现 ， 重 力 波 
的 幅度 并 不 都 是 随 高 度 的 增加 而 增 大 的 。 

几乎 所 有 的 雷达 观测 结果 都 报道 了 中 层 顶 区 域 
重力 波 的 垂直 波长 A 与 周期 7, 之 间 有 类 似 的 系统 性 
关系 A. 7% 。 在 Illinois (40"10'N,88"10'W) ,Gardner 
and Voelz 得 到 了 A.=0.247” 波长 与 周期 之 间 存 在 的 
相关 性 外 ;后 来 Collins 等 在 分 析 Wind/Temperature 钠 
雷达 数据 的 基础 上 ,得 到 了 A 和.=0.427%” 的 关系 中 ;在 
Sao José dos Campos(23°N ,46"W)]) , Yang 等 报道 了 A.= 
0.427%”， 的 关系 外 。 尽 管 各 个 地 方 都 得 到 了 相 类 似 的 
参数 关系 ,但 是 观测 到 的 Cp 值 却 变化 较 大 。 对 于 
一 个 给 定 的 周期 为 3h 的 常见 重力 波 ， 在 Ilinois 观 
测 到 的 垂直 波长 约 为 4.6 km (Gardner 的 数据 则 为 
6.9 km) ,在 Sio Jos& dos Campos 观测 到 的 垂直 波长 为 
3.9km， 在 北京 观测 到 的 却 是 3.54 km。 这 说 明 重 力 
波 活动 具有 明显 的 地 区 特性 ， 在 不 同 地 点 观测 结 
往往 不 尽 相同 。 

重力 波 的 线性 饱和 理论 认为 重力 波 垂直 波 数 谱 
的 形状 和 幅度 值 是 由 风 剪 切 和 对 流 不 稳定 性 控制 
的 , 当 重 力 波 的 幅度 增长 到 对 流 或 剪 切 不 稳定 点 时 ， 
重力 波 将 出 现 饱和 ， 饱 和 重力 波 的 垂直 波 数 谱 在 高 
波 数 区 域 的 谱 指 数 近似 为 n=-3 外 。Gardner 的 扩散 
滤波 理论 认为 波 致 涡流 扩散 是 重力 波 耗 散 的 主要 途 
径 ;在 重力 波 向 上 传播 时 会 被 背景 大 气 严 重 吸收 , 重 
力 波 垂直 波 数 谱 相 应 地 出 现 选择 性 滤波 效应 。 只 有 
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那些 有 效 垂直 
人 
中 层 顶 区 域 而 被 雷达 观测 到 ， 其 它 重 力 波 将 会 被 严 
重 衰减 而 消失 。 激 光 雷 达 在 中 层 顶 区 域 所 观测 到 的 
重力 波 都 是 接近 饱和 的 重力 波 (8 趋 近 于 0) ,幅度 也 
接近 最 大 值 ， 且 都 大 致 治 着 扩散 衰减 极限 线 w=Dki 


扩散 速度 (上 KD) 小 于 垂直 相 速 度 (w/) 
传播 到 


分 布 。 所 以 ,重力 波 参 数 之 间 示 具有 系统 相关 性 , 即 
A.x7T% KEcxc 和 KEc 记 20 在 北京 观测 到 的 重 
力 波 参 数 关 系 为 A.=0.2267T%”” KE=6.66x10-0K ”和 
KE=8.97x10“ 记 ,这 些 参 数 关 系 和 扩散 滤波 理论 、 
重力 波 的 线性 饱和 理论 所 预测 的 结果 是 一 致 的 。 
Gardner and Voelz 认为 波 致 涡流 扩散 是 重力 波 


耗 散 在 大 气 中 的 主要 途径 ， 涡 流 扩散 系数 D 可 以 通 
过 雷达 的 观测 数据 进行 计算 外 . 
D(1-2HB) FF A a 9) (5) 


其 中 ,4 为 由 于 粘 滞 效应 引起 的 衰减 系数 ;d, 为 传 
导 引 起 的 衰减 系数 .在 认为 波 致 涡流 扩散 效应 对 应 
的 衰减 系数 d, 远大 于 分 子 粘 湛 和 热传导 对 应 的 衰 


减 系数 ( 即 4,>>d,+ 人 )， 且 雷 达 在 中 层 顶 区 域 观测 


到 的 都 是 接近 饱和 a 趋 近 于 0) , 式 (5) 可 
以 简化 为 : 


” 0 (9) 

为 了 验证 Gardner 的 扩散 滤波 理论 ,图 6(a) 给 出 
了 北京 雷达 观测 到 的 162 个 重力 波 的 (Kk ,w) 统 计 分 
布 图 ;其 中 实 线 是 对 实验 数据 的 最 小 二 乘 拟 合 ,虚线 
是 由 式 (6) 计 算出 来 的 扩散 衰减 极限 线 w=Dk。, 图 中 
显示 绝 大 数 所 观测 到 的 数据 点 都 满足 w 三 Dk 的 条 
件 ， 即 雷达 观测 到 的 重力 波 都 满足 有 效 垂直 扩散 速 
度 小 于 重力 波 运 动 的 垂直 相 速 度 ， 是 没有 被 背景 大 

气 选 择 性 衰减 的 重力 波 。 

6(b) 是 雷达 观测 到 的 重力 波 的 幅度 值 (KE 动 
能 ) 随 w-Dk 有 的 分 布 图 ,图 6(c) 是 这 些 重 力 波 的 幅度 
增长 因子 6 随 w-Dk 有 的 分 布 图 。 幅 度 分 布 统计 结 
显示 ， 实 验 观 测 到 的 这 些 重力 波 并 不 是 完全 分 布 在 
w=D 甩 衰减 极限 线 附 近 是 整体 地 向 w>Dk 方 向 
偏 移 。 这 表明 雷 达观 测 到 的 重力 波 并 不 完全 是 幅度 
值 达 到 最 大 的 重力 波 ， 即 不 完全 是 饱和 了 或 趋 近 饱 
和 的 重力 波 。 
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(c) 重力 波幅 度 增长 因子 B 随 w-D 中 的 分 布 图 
(c) Distribution of amplitude growth factor of gravity waves versus w—DkR: 
图 6 北京 上 空中 层 顶 区 域 激 光 雷 达观 测 到 的 重力 波 参 数 统计 分 
析 图 


Fig.6 Statistic analysis on the parameter distribution of gravity waves 


observed by lidar in the mesopause region over Beijing 


图 6(c) 的 统计 结 一 步 显 示 , 在 满足 w 二 Dk? 
条 件 的 重力 波 中 ， 并 不 是 所 有 重力 波 的 幅度 都 是 随 
高 度 增 加 而 增 大 的 ， 还 有 部 分 重力 波 的 幅度 是 衰减 
的 (增长 因子 6 是 负 值 )。 图 中 国 圈 中 的 数据 点 明显 
远离 需 位 置 线 ( 即 w=D 尼 虚线 和 B=0 虚线 ) ,即使 雷 
达 数 据 统 计 分 析 中 存在 部 分 误差 (通常 小 于 3%), 这 
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些 数 据点 也 是 必然 大 量 存 在 的 。 这 说 明 Gardner 的 扩 
散 滤波 理论 中 仅 认 为 波 致 涡流 扩散 效应 是 重力 波幅 
度 衰减 的 主要 因素 是 不 正确 的 ;在 实际 大 气 重 力 波 传 
播 过 程 中 ,大 气 中 的 分 子 粘 滞 和 热传导 等 效应 也 能 

致 重力 波幅 度 的 衰减 ,它们 也 可 能 在 大 气 重力 波 的 饮 
和 和 耗 散 机 制 中 起 到 重要 作用 。 


4 结 论 
基于 北京 钠 激光 雷达 对 中 层 顶 区 域 重力 波 活动 
的 长 期 观测 ， 采 用 Yang 等 的 重力 波 人 参数 提取 方法 ， 


提取 了 162 个 钠 层 单 色 重 力 波 的 垂直 波长 .周期 . 幅 
度 值 和 幅度 增长 因子 ， 统 计 分 析 了 北京 上 空中 层 顶 
区 域 重 力 波 的 参数 关系 ,结果 显示 ,激光 雷达 所 观测 
到 的 重力 波 的 参数 具有 系统 相关 性 A.=0.2267W%™ 、 
KE=6.66x10 习 ki” 中 和 KE=8.97x10 六 及 ”“”。 这 些 参 数 
关系 和 Gardner 的 扩散 滤波 理论 预测 的 结果 是 一 致 
的 ,说明 由 于 重力 波 在 向 上 传播 的 过 程 中 ,背景 大 气 
对 重力 波 具 有 了 明显 的 选择 性 滤波 效应 ， 激 光 雷 达 在 
中 层 顶 区 域 观测 到 的 大 气 重力 波 大 致 满足 A.cc7%5、 
ee KEcxfi 参数 关系 。 但 是 ,统计 发 现 激 光 

达 在 中 层 顶 区 域 观测 到 的 重力 波 并 不 是 都 趋 近 饱 
0 志 生 本) 汪汪 市 站 苦力 下 的 十 庆 全 是 
减 的 ; 这 表明 扩散 滤波 理论 中 仅 认 为 扩散 效应 是 重 
力 波幅 度 衰 减 的 主要 因素 是 不 正确 的 。 
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